Edge States on a Quantum Hall Liquid— Solid Interface 
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We study the edge-states excitations of a droplet of quantum Hall liquid embedded in an electron 
(Wigner) solid. The presence of strong correlations between the liquid and solid sectors in the 
ground state is shown to be reflected in the density of states D(E), associated with the excitations 
of the liquid-solid interface. We find that the prominent effect of these correlations is a suppression 
of D(E) with respect to its value (Do(E)) in the absence of the Wigner solid environment: D(E) ~ 
e~ a ' E 'Do(E). The coefficient a (which is shown to vanish for a perfectly regular distribution of 
electron sites in the solid), is evaluated for two different realizations of an irregular distribution. We 
conclude that probing this effect (e.g. in a tunneling experiment), can provide evidence for correlated 
liquid-solid mixture states in quantum dots, or disordered samples, in very strong magnetic fields. 
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I. INTRODUCTION AND PRINCIPAL RESULTS 

The two-dimensional electron gas (2DEG) in strong 
perpendicular magnetic fields can form a variety of ex- 
otic quantum phases. In particular, in clean systems at 
moderately low filling fractions (close to v = 1/5), the 
correlations which favor a fractional quantum Hall liq- 
uid (QHL) state compete with the crystalline order of a 
Wigner solid (WS) Q^]. This competition can induce 
transitions between the QHL state and the insulator, as 
has been observed experimentally In the presence 
of slowly varrying disorder or a confining potential, the 
electronic ground state may develope a fractured order 
- namely, form a binary liquid-solid mixture. A liquid- 
solid separation possibly occurs also at higher filling frac- 
tions (in which the electrons in the solid regions should 
form a glassy state controlled by short-range disorder). 
An experimental evidence for this scenario is provided by 
photoluminescence data j|. In addition, a set of puzzling 
transport data |^|| can be explained most naturally un- 
der the assumption of a macroscopic inhomogeneity 0. 

The interplay of QHL correlations and crystallization 
in the low v regime has been clearly demonstrated by 
Zheng and Fertig 0j. Using a variational calculation, 
they have shown that a Wigner lattice with an interstitial 
electron introduced via a Laughlin-like Jastrow factor, 
can be lower in energy compared to the perfect WS with 
the same total number of electrons. This implies that 
in a certain range of z/s, the crystal is unstable to a 
specific type of density fluctuations - pre-formed QHL 
droplets. In the presence of density fluctuations induced 
by a slowly varrying external potential, it is reasonable 
to expect a nucleation of such interstitials in the higher 
density regions. It is therefore suggestive that the ground 
state slightly below v — 1/5 separates into QHL and 
WS sectors, which are correlated by a Jastrow factor to 
minimize the energy of electrons close to the liquid-solid 
interface (see Eq. (|])). 



In the present paper we investigate the physical im- 
plications of liquid-solid mixture states, as reflected by 
the corresponding low lying excitations. Similarly to a 
finite droplet of a primary QHL (of v = 1 /m with m an 
odd integer) , the gapless excitations are chiral edge states 
Q, i.e., deformations of the boundary of the incompress- 
ible droplet, which travel in a definite direction along 
the boundary. However, in case the liquid droplet is em- 
bedded in an electron solid (rather than a vacuum), the 
nature of these excitations of the liquid-solid interface is 
affected by the correlations between the two sectors. In 
particular, high amplitude deformations of the interface 
are generally suppressed, since the liquid electrons are 
constrained by their tendency to avoid the proximity to 
localized sites of the WS as much as possible. This can 
lead to a decay of the density of states with increasing 
deviation from the Fermi level, as long as higher energy 
excitations (which involve, e.g., a reorganization of the 
WS electrons) are not yet activated. 

To facilitate the derivation of this peculiar effect, we 
consider a simple geometry of a large, circular quantum 
dot, in which the electrons are assumed to form a disc 
of QHL surrounded by a WS (see Fig. 1). We evaluate 
the electron propagator and consequently the density of 
states for tunneling into the liquid-solid interface, D(E). 
In the thermodynamic limit, we find 



D{E) ~ e~ alEl D {E) 



(1) 



where Do(E) corresponds to the ordinary edge-state (on 
an interface between QHL and a vacuum) , and the coeffi- 
cient a depends crucially on the distribution of localized 
sites in the WS sector. In particular, when these sites 
form a structure with a perfect crystaline order arround 
the disc (and the lattice constant is commensurate with 
the circumference), a = 0: in that case, the WS electrons 
merely deform the effective boundary of the liquid into 
a regular shape, as depicted in Fig. 1(a). In contrast, 
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an irregular distribution of sites induces frustration, and 
thus a suppression of D(E). 

The most direct way to probe this effect is via tunnel- 
ing into the 2DEG, e.g. using the technique developed 
by Ashoori 0]: the tunneling conductance is given by 
G(V) ~ D(eV) (where V is the voltage across the tun- 
nel barrier). The suppression of G(V) at low V may 
lead to a non-monotonous behavior - at higher voltage 
bias, higher energy excitations take over and induce an 
increase of G(V). The effect is expected to become more 
pronounced with increasing inhomogeneity of the exter- 
nal potential. In particular, in a quantum dot where 
both the number of electrons and the confining potential 
can be controlled, the suppression of G(V) is expected to 
exhibit oscillations: the tunneling rate should be maxi- 
mized when the control parameter enables a nearly reg- 
ular configuration of sites in the WS sector. 

In the following sections, we detail the derivation of the 
electron propagator along the liquid-solid interface (Sec. 
, and the implied behavior of the density of states, Eq. 
(Sec. III). In the latter, we consider two different 
realizations of the irregularity in site configuration: (a) 
a regular distortion of the circular symmetry, and (b) a 
symmetric random distribution. The corresponding ex- 
pressions for the su ppr ession-time a are given by Eqs. 
(If) (case (a)) and @,@ (case (b)). 



II. DERIVATION OF THE ELECTRON 
PROPAGATOR 

As we explained in the introduction we expect that 
the ground state wave function that describes the QH 
droplet surrounded by WS electrons is essentially of the 
following form: 

*ls = A{y L (zi, ■ ■ ■ , z N )V s (wi, . . .,w M ) Y[( Zl - Wj) m } 

(2) 

where ^ l describes the liquid part, ^ s the solid part, and 
the last expression describes the Jastrow correlations of 
these two phases. A denotes antisymmetrization over all 
electron coordinates, m corresponds to the filling factor 
v = 1/m of the QH liquid part: 

..., ZN ) = H( Zl - Zj r ex P {-i J2 n 2 } ( 3 ) 

i<j 

The wave function 
9 s (wi,...,w M ) =exp{i^uJ l i? i - ^J2\ m \ 2 } ( 4 ) 

describes WS electrons localized on the positions Ri,i = 
1, . . . , M. (It is a multiple of lowest-Landau-level delta 
functions.) 

Our main assumption in the derivation of the electron 
propagator is that the edge of the QH droplet, coupled to 



WS, essentially behaves as a slightly modified Luttinger 
liquid. That assumption allows us to use a construction 
of one electron state on the edge similar to the one of 
the Luttinger liquid. At the end of the derivation we will 
be able to specify constraints on the configuration of the 
WS sites, such that the assumption is valid. 

Because of the above assumption, it would be instruc- 
tive first to briefly recapitulate the derivation of the 
(equal-time) propagator when the QH disc is surrounded 
by vacuum || . 

The derivation begins by considering m - Laughlin 
quasihole constructions, i.e. 



A' 



(5) 



where £ lies outside of the system; |£| > R, and R = 2mN 
is the radius of the QH droplet. \l/z, is the Laughlin wave 
function (Eq. (||)). The first step towards the electron 
correlator is calculating the following scalar product of 
the state (|) 



-2mN 



(6) 



The leading contribution in the R[\£\ expansion can be 
obtained by the plasma analogy H or simply by consid- 
ering the decomposition of coordinates into the center of 



mass, 



Z C m — jf Yli=i z ii an d relative ones. When the 



numerator is approximated by 



>]Tln|z fe -£| 2 



Nm2 In |£| - mN^ - mN 



z* 



(7) 



the integration over the the center of mass coordinate 
is decoupled from the other integrations which do not 
depend on £. It yields the leading dependence on -R/|£l : 



(8) 



In an average, macroscopic picture, we expect singular 
behavior as |£| — * R, with singularity at |£| = R (for the 
equal-time, equal-space correlator) ||. Then Eq. (j|) can 
be rewritten as 



^(r,o = (i-|j5) 



(9) 



and we will assume that it is valid also for |£| ~ R (|£| > 
R). By doing this we neglect any finite size corrections 
(due to finite R) which might be present in J\f a as |£| — > R. 
To get the electron propagator we analytically continue 
the function Af Q (£* , £) of the variable £ to J\f (£* , £) which 
depends on £ and £. That allows us to take £ and £ to 
the edge of the system - £ = i?exp{z2-7r^} and £ = R, 
without encountering the singularity for x = (which 
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determines the behavior of the function in its neighbor- 
hood) . By taking £ and £ to the edge, we in fact describe 
a particle-hole excitation on the edge that goes into the 
electron propagator, and find that for x <C L the propa- 
gator behaves as 



g°(x) = M,(r,or (JV-1)m £ (JV-1)m 

Jxyl exp{«m(iV - 2)-^} 



(10) 



m(N — |)^r in the exponential is the value of the gener- 
alized (m 7^ 1) "Fermi momentum". 

Following a similar strategy, we address the case of 
the droplet surrounded by WS, described by the wave 
function (0). We consider the following state 



N 



*(0 = Hi** - Q m *L S 



and the corresponding scalar product: 



(11) 



(12) 



where 



y ]^[d 2 Zj y ]~[ c? 2 Wj exp{y~^ 2to ln|zj 

»'=i j=i *<i 

n a; ^ n 

+EE 2mln i z, "^i ~ 2 E ' Zi ' 2 

i=l j=l i=l 
M JV 

-,Z)K-i2j| 2 + 2mX;ia|« i -e|} 



(13) 



3=1 



In the above formulas the antisymmetrization was ne- 
glected, which is possible due to the localized nature of 
the WS electrons. Again, if for If | ^> R and \wj\ ^> 
R, j — 1, . . . , M, the last sum in Q13) is approximated as 
(0), and 



AT M 



^^2mln|zi- wj\ « 

i=i j=i 

M M M 

J2 2mN In \ Wj \ - mNZ cm ^ mNZ* m £ — , (14) 

^ ' ^ ■* -7/ ■ * 111 



. wz 
3 3 



the integration over Z cm yields the functional dependence 
of the integrals over z^s in (]l^) on £ to leading order in 
ft I Sj=i S~ + 1 1- It is of the following form: 



M 



d + ^ 2 iE^ + 7i 2 )- 



(15) 



Now we will assume that the sites Rj , j = 1, . . . , M 
are such that the integration over Wj's is dominated by 
contributions for which 



M 



ft 2n iE^r + -i*- 



n = 1, . 



(16) 



Then, the result of the z-integration can be expressed 
as an expansion in variables symmetric in £ and Wj's, 
defined as 



M 

E 



1 1 



(17) 



The first two terms of the expansion are given by the 
expression (15). As explained in Appendix A, the final 
w-integration amounts to replacing Wj with Rj, defined 



in the Appendix, and variables (17) in the expansion yield 



M 

X n {R j ,j = l,...,M}=J2-&,, 



R n 

3=1 n 3 



1 

¥ 



(18) 



At this point, we can see that in the case where the po- 
sitions of WS electrons satisfy 



M M 

E — =E 



1 

R* n 



0, n= 1, 



, 00 



(19) 



the problem reduces to the one of the droplet surrounded 
by vacuum, and the expansion should sum up to the 
Luttinger- liquid form (||) . These conditions can be satis- 
fied, e.g., when 



Rj = Rexp{i(Qj +d a )} 



(20) 



where 9j = j9, 9 = -A, and R is a constant radius, i.e., 
when a commensurate chain of WS electrons surrounds 
the droplet (see Fig. 1(a)). Note that our assumption 
of a small correction to the Luttinger-liquid behavior is 
justified, provided the configuration of WS sites is a small 
perturbation of one which satisfies Eq. (|l^) . 

To get the electron propagator, we first exponentiate 
the expression of the leading order behavior for |£| ^> 
R, \Rj\^>R, j = l,...,M 



M 

AT(r,o«(i-H 2 iE 5- 

1 R-i 
3=1 3 



(21) 



Here we assume that even for Rj's close to the droplet, 
Y^, -£7 is small with respect to t|t ~ , (which implies 

that the pole structure of the correlator is similar to the 
one of the Luttinger liquid with the pole slightly shifted 
from |£| = R). Hence, we can also regard Eq. (|2l| ) valid 
for general Rj, j = 1, ...,M, and £ (with |i?j|,|£| > 
R). A/"(£*)£) is then analytically continued to A/"(£*,£); 
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taking £ = i?exp{i27r ^} and £, = R gives for the electron 
equal-time propagator 

G e (x) ~ exp{im(N~ ±)27rf } X 
[exp{i7r^} — exp{— i7r^} 
- exp{i7r^}£ - exp{-i7r^}S* 
-exp{i7rf}|E| 2 ]-" 1 



G e (x = 0, t > 0) = -i < 0|*(0, t)* t (0, 0)|0 >, (26) 



(22) 



where £ = ££i # 



In the limit £ <C 1 that corresponds to the short distance 
behavior, the electron propagator can be expressed as 



1 x 1 
G e {x) ~ exp{im(A - -)2tt-} x — — 

2 L (is + C\x — C-2) 



where 



and 



Ci 



2Z{£} 



C 2 = - 



L(2ft{£} + |£| 2 ) 



(23) 



(24) 



(25) 



(TZ{T,} and !{£} are the real and imaginary parts of £, 
respectively). Note that |£| serves as a small parameter 
that determines the deviation from a standard Luttinger 
liquid behavior (Eq. ([To|)). 



III. THE TUNNELING DENSITY OF STATES 

In the previous section we calculated the equal time 
propagator. In order to get the tunneling density of 
states, we need the time-dependent propagator G e (x = 
0, t). In general, this requires the knowledge of the ener- 
gies of excited states of the system. For WS configura- 
tions that satisfy the condition s ( f l9|) (and hence have the 
Luttinger liquid correlations (|l0|)), we expect that the 
ground state is one of the edge states of a QH droplet 
surrounded by vacuum; the excited states (which may 
also be interpreted as edge states of that droplet) have 
energies linear in momentum, measured from the new 
ground state. The linear dispersion is expected on gen- 
eral grounds, as the first order expansion in small mo- 
menta, and not precluded by any (symmetry etc.) argu- 
ment fiTo] . Then, to get G e (x — 0,i) we should merely 
substitute the coordinate x with vt (where v is the veloc- 
ity of the drift motion of electrons on the boundary || ) . 
The sign of time should be specified, as we will explain 
and elaborate below. Assuming that the linear dispersion 
is valid also in the case of a small deviation from Eq. jl9| ) 
(i.e.. for |S| <C 1), we may apply the same substitution 
in (§3|) to get G e (x = 0,t). 

We recall the definition of the (equal-space) fermionic 
Green's function, in the ficld-thcoretic notation: 



and 

G e (x = 0,t < 0) = i < O^O, 0)*(0,f)|0 > (27) 

where VE^ and ^> are the electron and hole creation op- 
erators, respectively. If we assume that, in our system, 
equation 

< 0|*(0, i)* f (0, 0)|0 >=< 01^(0, t)tf (0, 0)|0 > (28) 

holds, as it does in the case of the standard Luttinger 
liquid, going from formula ( p7[ ) to formula (26) involves 
only time translation (T — > T — t) and time inver- 
sion (T — > — T) (and the overall change of the sign). 
In our case, the particle coordinates £ = Rexp{ix^-} 

and £ = R become under the substitution (x — > vt) 
£(t) = Rexp{ivt^-} and |(0) = R. £(i) and £(0) de- 
note the hole and electron coordinate respectively. If we 
reconsider the correlator £) with the substitutions, 

we will get for the hole propagator ( p7|) in the short time 
limit: 



G e {0,t < 0) oc 



1 



(ivt + Civt - c 2 y 



(29) 



On the other hand, to get the electron propagator for 
t > (|2^) we should perform the time translation and 
inversion. This amounts to an exchange of £(t) and £(0), 
which leads to the following short time behavior: 



G e (0,t > 0) cx 



1 



(-ivt + Civt - c 2 y 



(30) 



If we consider t as a complex variable and assume |S| <C 
1, we may approximate the positions of poles in (E9) and 
( |30| ) as ti ~ —i^ 2 -, and t% ps i— , respectively. Then, the 
tunneling density of states D(E) is given by 

D(E)-n { j\teMiEt} i . t _ 1 C2/v)m + 

f°° l 

/ dt exp{iEt}- ; - } (31) 

Jo J (it + C 2 /v) mf v ' 

where E is the energy measured from the Fermi energy, 
ft = 1. As a final result we get 

IC I 

D(E) oc exp{-a\E\}D a (E) , a = — , (32) 

v 

where D (E) is the Luttinger liquid tunneling density 
of states. Therefore, to lowest order in |£|, the domi- 
nant modification to the standard Luttinger behavior is 
the exponential suppression, at a time scale a of order 
~ (i/u)|£|. It should be stressed that a, and hence 
D(E), is a local quantity (adiabatically varrying around 
the disc) that describes the tunneling density of states of 
an electron at distance ~ R from the origin, and at angle 
7T-j <C 1 to the reference point. 
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Below we calculate a for the two different types of 
imperfect site configurations depicted in Fig. 1(b), (c). 

A. Inhomogeneous Configuration of Sites 

We model the inhomogeneous configuration (Fig. 1(b)) 
by considering M — 2n + 2 WS electrons, at a distance 
R from the origin, where two of them are exactly on the 
opposite sides of the droplet, i.e. the sum of their phases 
is exp{i7r} + 1 = 0. The rest 2n electrons are positioned 
at the angles 9j and —9j (1 < j < n), where 



9j = 0j + ecos{S{j-l)} , 



(33) 



9 = ^pj- and S = e > represents a small deviation 
from the perfect chain distribution, which is modulated, 
as we move from the reference point at angle 7 = 0. We 
then get S = S , where 

R n 

S ° = $Z ex P{^i + «ecos{(5(j - 1)}} + 
R ■ , 

3=1 

R n 

^ V exp{-idj - ie cos{5(j - 1)}} . (34) 
R i=i 

For e small So can be approximated as 
R n 

S ~ V exp{i0j'} cos{S(j - 1)} 

R . , 
j=i 

R " 

_ ie V exp{-i0j} co S {5(j - 1)} (35) 
« i=1 



and for n ^> 1 this yields 



8nei? 
°~~3^S 



(36) 



If the reference point is at an arbitrary angle 7^0, 
the sum becomes 



phase of the Rj's; here we discuss the effect of each type 
of randomness separately. 

In a phase-distorted chain at radius R, the site Rj is 
given by 



Rj = Rexp{i(j6 + Sj)} 



(39) 



where 9 — and Sj is a random variable. Averaging 
over the distribution of <5j's, we obtain 

M 

<£>=£<f> = 



A/ 



|<exp{-z<5})^exp{-i(j0)} = O 
R 3=1 



(40) 



The lowest order contribution to Ci (Eq. (|2q)) is then 
- (|S| 2 ), where 



s i 2 ) = £E( 2 - 2 ( ex pM))- 



i? 2 



(41) 



For a Gaussian distribution P(<5) = (l/v27rer)e^ 
we get 

E>2 r>2 Mrr 2 

<|S| 2 ) = - exp{-a 2 /2}) « |^ , 



(42) 



where the last approximation holds for u < 1. Substi- 



tuting in Eqs. ([25 



and (|32|), this yields 

L R 2 Ma 2 
~ 2tto^2 2 ' 



(43) 



We next consider a radius-distortion of Eq. (pp| ) of the 
form 



i?j = (fl + r^exp^)} 



(44) 



E = S cxp{-«7} . 



Inserting in Eqs. (E5h and 



we then get 



L 

7TV 



2S cos{ 7 } + |S | 



(37) 



(38) 



Therefore, at 7 = 5, a becomes a oc e 2 , at 7 = it 
is a oc (e — Const e 2 ), Const > 0, and, at 7 = tt, it is 
a oc (e + Const e 2 ). This means that around 7 = and 
7 = 7r deviations from the perfect chain case are stronger, 
and the suppression of the density of states is stronger in 
the denser region (7 = tt), but only to second order in e. 



B. Random Distribution of Sites 

We consider a configuration of sites obtained by a ran- 
dom distortion of the perfect chain Eq. ( pp| ) (see Fig. 
1(c)). A general distortion modifies both the radius and 



where rj is a random variable, subject to a distribution 
of width oy -C R. Again, the first order in S vanishes 
upon averaging: 



R 3-1 



(45) 



For a symmetric distribution of r^'s, (r) = and we get 
?2 M p2 



(46) 



We then get an expression for a which is quite similar to 
Eq. <M): 



L R^ 1 o 
a re — -p-Mcr 2 



(47) 



5 



ACKNOWLEDGMENTS 

We thank E. Akkermans, R. Ashoori, Y. Avron, J. 
Feinberg, Y. Meir and S. Sondhi for fruitful discussions. 
This work was partly supported by grant no. 96-00294 
from the United States-Israel Binational Science Founda- 
tion (BSF), Jerusalem, Israel, and the Technion - Haifa 
University Collaborative Research Foundation. M.M. 
also acknowledges support from the Fund for Promotion 
of Research at Technion, and the Israeli Academy of Sci- 
ences. 



APPENDIX A: 

To evaluate < > (Eq. @), one needs to 

solve integrals of the form 



<fiwf(w,w*)\w\ 2mN e-i lwl2 ei wR hi w * R > , (Al) 



where w is the coordinate of the j'th WS electron (the 
index j being omitted), and f(w,w*) can be expressed 
as a power series in R/w, R/w*: 



f(w,w*) = V* a nk ( — 

n,k=0 v 



(A2) 



Eq. QA1[ ) is then recast as 

oo 

ImN —n,mN —k ) 



n,k=0 



I p .i = I d 2 wvfw* l e-i\ w \ 2 ei wR iei w * R i 



(A3) 



Below we show that provided (mN — n), (mN — k) ^> 1, 

R 2 



ImN—n,mN—k ~ IrnN ,mNRj " (Rj , 



= -j = e a ^ R ^ where a(\Rj\) 
Rj R* 



\R 3 \ 2 



(A4) 



note that exp{a(|i?.,|)} — > 1 for (R/Rj) <C 1. In addition, 
we note that since the integration over w is domi nate d by 
w ^ Rj, and Rj > R, the series expansion Eq. (A2) can 



be cut at some n c , k c such that \R/Rj\ n ", \R/Rj\ kc <C 1, 
yet (mN — n c ), (mN — k c ) 3> 1. Consequently, we obtain 

I ~ f(Rj,Rj)ImN,mN , (A5) 

which implies that the {w ? }-integrations over the terms 
X n (Eq. @)) yield Eq. @. 

We now derive the approximation Eq. ( |A4| ) - the cen- 
tral result of this Appendix. The integrals I P} i (Eq. (|A5|)) 
can be expressed as 



dR f dR\ 



wR i e i w " R j 



= 2 P+ 



8R? dR 



[WH^ 1 '] . (A6) 



A straight-forward application of the derivatives then 
yields 



J p ,,=47re4l^l R P R* l S(p,l) 



min{p,/} 



where S(p, I) = /J 



p\l\ 



^ i\(p-i)\(l-i)\ \\Rj\ 



(A7) 



We next assume that the above sum is dominated by 
1 <C i -C p, I, so that the factorials are well approximated 
by Stirling's formula: 



pie 2 



/ /?,I J V 1 p) V 1 



(A8) 



pi 



V 2m(p — i)(l — i) 
The sum is then replaced by an integral 



S(p,l) 



1 



(x) = x < In 



2ple 



2tt 

;J('-t) 

+0(ln(a;)) 



(A9) 



which can be solved in a saddle-point approximation. 
The saddle-point equation 4>'(x) = implies 



2ple 

Jr~S 



'-DC 1 -?)--*! 1 



p — X I — x 

(A10) 



and hence 



exp < x s 1 



(All) 



where x s is the solution of Eq. (ALC). For (R/Rj) <C 1 
(where p, I < Nm = R 2 /2) 



\R. 

and we get 

S(p, I) « ^/x~ s e Xs 



2 P l , jjii \ 1 

— s «p,J, 0"(x s )« — 



••' .S' 



2p2 



exp 



2p/ 



(A12) 



(A13) 



G 



Noting that p,l ~ R 2 /2, Eq. ( |A13[ ) implies Eq. (|A^). 
Q.E.D. 
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FIG. 1. A disc of QH liquid surrounded by a WS (the 
dots denote the sites at which the WS electrons are localized) . 
Three different types of site configurations are sketched: (a) 
an ordered, commensurate chain (the effective boundary of 
the liquid in the presence of the WS electrons is marked by 
a dashed line); (b) an inhomogeneous configuration; (c) a 
random distribution. 
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